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ABSTRACT 

The non-local thermodynamic equilibrium (NLTE) line formation of neutral 
boron in the atmospheres of cool stars are investigated. Our results confirm that 
NLTE effects for the B I resonance lines, which are due to a combination of 
overionization and optical pumping effects, are most important for hot, metal- 
poor, and low-gravity stars; however, the amplitude of departures from LTE 
found by this work are smaller than that of previous studies. In addition, our 
calculation shows that the line formation of B I will get closer to LTE if the 
strength of collisions with neutral hydrogen increases, which is contrary to the 
result of previous studies. The NLTE line formation results are applied to the 
determination of boron abundances for a sample of 16 metal-poor stars with the 
method of spectrum synthesis of the B I 2497 A resonance lines using the archived 
HST/ GHRS spectra. Beryllium and oxygen abundances are also determined for 
these stars with the published equivalent widths of the Be II 3131 A resonance 
and O I 7774 A triplet lines, respectively. The abundances of the nine stars 
which are not depleted in Be or B show that, no matter the strength of collisions 
with neutral hydrogen may be, both Be and B increase with O quasi-linearly in 
the logarithmic plane, which confirms the conclusions that Be and B are mainly 
produced by primary process in the early Galaxy. The most noteworthy result 
of this work is that B increases with Fe or O at a very similar speed as, or a bit 
faster than Be does, which is in accord with the theoretical models. The B/Be 
ratios remain almost constant over the metallicity range investigated here. Our 
average B/Be ratio falls in the interval [13±4, 17±4], which is consistent with the 
predictions of spallation process. The contribution of B from the //-process may 
be required if the 11 B/ 10 B isotopic ratios in metal-poor stars are the same as the 
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meteoric value. An accurate measurement of the 11 B/ 10 B ratios in metal-poor 
stars is crucial to understanding the production history of boron. 

Subject headings: Galaxy: evolution — line: formation — line: profiles — stars: 
abundances — stars: atmospheres — stars: Population II 



Introduction 



The origin and evolution of boron (and beryllium) are of special interest because they are 
hardly produced by the standard big ba ng nucleosynt h esis ( BBN), nor can they be produced 
by nuclear fusions in stellar interiors. iReeves et al.l (119701 ) were the first to propose that 
most of the light elements (Li, Be, and B) can be produced by the high energy processes 
involving Galactic cosmic rays (GCR) act ing on the interstellar medium (ISM). Subsequent 
detailed model by iMeneguzzi et al.l (ll97l[ ) could reproduce most of the observations at that 
time; one of the exceptions is that the theoretical isotopic abundance ratio 11 B/ 1O B0 (2.4) is 
lower than the observed solar meteor itic value of about 4.0. In order to solv e this problem. 



they i ntroduced a non-observable intense low energy component to the GCR. IWoosley et al. 
(Il990h proposed that 7 Li and 11 B can be synthesized by neutrino spallation in the helium and 



carbon shells of supernovae, i.e., the so-called //-process, which argues against the existence of 
an unobserved low energy component of GCR. Nevertheless, the GCR process, which predicts 
a quadratic relation between B (also Be) and O abundances, remained as the standard 



(Duncan et al. 


1992; 


Edvardsson et al. 


1994; 


Duncan et al. 


1997. 


1998b: 


Garcia Lopez et al. 



19981 ). as well as on Be in metal-poor stars showed that B and Be abundances increase 
quasi-linearly with O abundances, which indicates that B and Be are probably produced by 
primary process, rather than the standard secondary GCR process. 

All the studies on B a bundances menti oned above were based on the ultraviolet (UV) B I 
resonance lines at 2497 A. iKiselmanl ( 19941 hereafter K94) investigated the formation of the 
B I resonance lines in the atmospheres of solar-type stars and found significant departures 
from LTE in stars ho t ter or more metal-deficient than the Sun. In their subsequent work, 
Kiselman fc Carlsson (jl996 . hereafter KC96) presented the NLTE abundance corrections 
for the B I 2497 A resonance lines for a grid of cool stellar model atmospheres. After 
that, in almost all the studies on boron, B abundances were determined in LTE first and 
then corrected with the results of KC96. However, NLTE corrections are model dependent 



1 A/B = N(A)/N(B) 
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and in pri nciple should be a pplied only to the analysis using the same model atmospheres. 



Moreover, iBoesgaard et al.l (120041 ) pointed that B abundance with the NLTE corrections of 
KC96 increases more slowly for halo stars than the Be abundance does, which is not predicted 
by light-element synthesis or depletion, therefore they suggested a full NLTE analysis for B, 
rather than simply applying corrections to the LTE abundances. 

The aim of this work is thus to re-investigate the NLTE line formation of neutral boron 
in the atmospheres of cool stars first, and then to obtain the NLTE B abundances for a 
sample of metal-poor stars, based on which the origin and evolution of B in the early Galaxy 
can be investigated. In § 2 we briefly describe the observations and data reduction. Section [3] 
presents the adopted model atmospheres and stellar parameters. Section 0] deals with the 
NLTE line formation of neutral boron. Section [5] gives the abundances and uncertainties. In 
§ If)] we discuss the implications for the origin and evolution of boron, while in the last section 
we briefly summarize our results and conclusions. 



Observations and data reduction 



Our analysis of B abundances are based on the archived spectra of 16 cool metal- 
poor stars (—2.7 < [Fe/H|§ < —0.5), which were obtained with the Goddard High Resolu- 
tion Spectrograph (GHRS) on board the Hubble Space Telescop e (HST). The sampl e stars 
were selected from four original studies on boron abundances ( iDuncan et al.l Il992l Il997l ; 
Garcia Lopez et al.lll998t iPrimas et al.lll999l ). All the stars were observed using the G270M 
grating centered at the B I 2497 A resonance doublet with a resolution of about 25,000. The 
signal-to- noise ratios (S/N) of the spectra are higher than 50 per diode. More details about 
the observations can be retrieved from the HST Science Archive Database and the afore- 
mentioned original paper. The spectra were reduced following the standard HST procedure 
using the STSDAS package of the IRAF0 suite of programs. 



2 [A/B] = log[JV(A)/J\T(B)]* - log[tf(A)/JV(B)] 

3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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3. Model atmospheres and stellar parameters 



We adopte d the opacity sampling (OS) model atmosphere MAFAGS-0S9 presented 
by iGrupp et al.l (120091 ) both in the NLTE line formation investigation and in the abundance 
analysis. This model is base d on the one-d imensional (ID) LTE t heoretica l mode l atmosphere 
MAFAGS-OS described by IGrupp! (120041 ) but incorporates the iKuruczl (120091 ) iron atomic 
data as well as the new bound-free (&-/) opacities for Mg I and Al I. The emergent solar 
flux from this new model shows significantly better agreement with observations, es pecially 
in the UV region of the spectrum, as can be seen in Figure 3 of IGrupp et al.l (120091 ) . 



Stellar parameters are fundamental to deriving elemental abundances in the photo- 
spheres of stars. However, it's impractical to determine them accurately with the narrow- 
wavelength-range (about 40 A) and severely-blended UV HST/GBRS spec tr a. We note d 



that 1 1 of our sample stars had b een s tudied prev i ously by lKorn et al.l (120031 ). iGehren et al. 



( 120061 ) , iMashonkina et al.l (120081 ) , and iTan et al.l (120091 ) , where the stellar parameters were 
determined spectroscopically using almost the same model atmosphere code and technique, 
i.e., effective temperatures were derived by fitting the wings of Balmer lines; surface gravities 
were based on the Hipparcos parallaxes ( Perryman et al. 1997 ) if availabl^l; iron abundances 
were determined from Fe II lines; and the microturbulent velocities were estimated by re- 
quiring the iron abundances independent of equivalent widths or by fitting the strongest 
Fe II, Ca I, and Mg I lines which are sensitive to microturbulent velocities. For the rest five 
stars, fortunately, high-resolution and high S/N spe ctra covering a wid e wavelength range 
were available from th e archived ESO VLT/UVES dDekker et aJboOoh or OHP ELODIE 
( IMoultaka et al.l 120041 ) spectra database. In order to keep our analysis consistent, we deter- 
mined the stellar parameters of these five stars using th e method describ ed above. The final 
stellar parameters are given in Table [TJ As discussed in ITan et al.l ( 120091 ) , the typical uncer- 
tainties for the effective temperature, surface gravity, iron abundance, and microturbulent 
velocity are estimated to be ±80 K, ±0.15 dex, ±0.08 dex, and ±0.2 km s -1 respectively. 



4. NLTE line formation of neutral boron 

4.1. Boron atomic model 

The boron atomic model adopted in this work is illustrated by its Grotrian diagram as 
shown in Figure [TJ Details for each item, as well as a comparison with the counterpart of 



4 Only one star BD — 13°3442 has no Hipp arcos parallax data; its su rface gravity was determined from the 
ionization balance between Cai and Can by Mashonkina et al. ( 20081 ). 
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K94 are described in the following paragraphs. 



4-1.1. Energy levels 



The energy levels were taken from the laboratory data compiled by lKramida &: Ryabtsev 



(120071 ) . Fine structure splitting was neglected except for the ground state 2p 2 P°. The final 
model atom includes 33 B I bound levels plus the B II ground state as the continuum. In 
fact, for testing purpose, we have extended our model by adding 61 B II bound levels and 
the B III ground state as the continuum, but the calculated results showed that there are no 
essential differences. The final boron atomic model adopted by K94 includes 30 bound levels 
and one continuum level, which is comparable to th at of this work. Th e energ y levels of K94 
were taken from the laboratory data compilation of iBashkin fc Stonerl (119781) except for the 



higher P° levels which were taken from the TOPbase (jCunto et al.lll993[ ). The mixing of 



laboratory and theoretical data makes the energy level sequence of K94 a bit different from 
that of this work, for example, their 7s 2 S level is higher than 7p 2 P°, a n d the 8s 2 S level is 
higher than 8p 2 P°, which are contrary to that of Kramida &: Ryabtsevl ( 120071 ). 



4-1.2. Radiative transitions 

Oscillator strengths were extracted from the TOPbase except fo r the transitions 2 p 2 P° - 
3s 2 S and 2p 2 P°-2p 2 2 D, which have accurate laboratory data from I Johansson et al.l ( 119931 ). 



K94 adopted exactly the same data as ours except that they excluded the lines with / < 0.001 
and the lines of extremely long wavelength. 

Photoionization cross-sections were available from the TOPbase for 26 levels; for the rest 
levels (the 4 P and 2 F° levels) hydrogenic approximations were adopted. K94 also adopted 
the photoionization cross-sections from the TOPbase when available, but their data for the 
4 P level was provided by Dr. K. Berrington and the 2 F° levels' cross-sections were copied 
directly from the 2 P° levels with similar energies. 



4-1.3. Collisional transitions 



For electron impact excitation, the effective col lision strengths were a vailable from the R- 
matrix with pseudostates (RMPS) calculations of iBallance et al.l ( 120071 ) for the transitions 
among the first 12 levels; the remaining allowed and forbidden tr ansitions were approxi- 
mated by the formulae of Ivan Regemorterl (119621 ) and lAUenl (119731 ) , respectively. Electron 
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impact i oniza tion cross-sections were calculated from the classical path approximation of 
Seatonl (119621) . K94 a d opted the electron impact excitation cross-sections calculated by 



Nakazaki fc Berringtonl (ll99ll) fo r the transitions among the first nine levels; the approxi- 



mation of Ivan Regemorterl (119621 ) was adopted for the remaining transitions. The el ectron 
i mpac t ionization rates of K94 were calculated from the hydrogenic approximations of I Allen 
jl973h . 

Excitation and ionization induced by i nelastic c ollisions wi th neutral hydrogen atoms 



are usually calculated using t he formula of iDrawinl (119681. 19691 ) presented in the form for 
ast rophysical app l icatio ns by ISteenbock _fc Holwegerl (1984 ). How ever, it has been shown 



by iBelyaev et al.l ( 119991 . for Na I) and iBelyaev &: Barkleml (120031 . for Li I) that collision 
rates of the resonance transitions given by Drawin's formula are overestimated by orders 
of magnitude. Therefore the collision efficiency was commonly constrained by applying a 
scaling factor Sh to the Drawin's formula empirically, and the factor cou ld be determined 



by requiring the deriy ed abundances from different lines to be consistent (IZhao et al.l 11998 



Zhao fc Gehrenl |2000| ). In this work, the cross-sections are computed with three different 
scaling factors, namely Sh = (no collisions with neutral hydrogen), 0.1, and 1, but as there 
are only two observed B I lines, and one is severely blended, we are not able to determine 
the final Sh value using the method described above. However, as will be discussed in § [61 
the uncertainty in the strength of collisions with neutral hydrogen does not affect our final 
conclusions much. For K94, the cross-sections for the b-f transitions induce d by collisions 
with n eutral hydrogen were calculated according to the analytic expressions of iKunc fc Soon 
(Il99ll ). while the bound-bound (b-b) transitions were neglected. 



4.2. NLTE line formation results 



The coupled radiative transfer and statistical equilibrium equations fo r boron were 
solved with a revised version of the DETAIL program (IButler fc Giddingslll9851) based on the 
accelerated lambda iteration following the efficient method described by lRybicki &: Hummer 
(Il99ll . Il992l ). In addition to the continuous background opacities, line opacities introduced 
by hydrogen and metal lines calculated with OS method were also taken into account in 
solving the equa tions. Lines shortw ard and longward of 2500 A wer e extracted from the 
VALD database (iKupka et al.lll999l ) and the line list of iKuruczl (119921 ). respectively. Please 
refer to IShi et al.l (120081 ) for more details about the dealing with background opacities. 



The departures from LTE for the populations of the first eight energy levels of B I 
and the ground state of B II calculated with different collision strengths with neutral hy- 
drogen are shown in Figure [2] for the representative metal-poor star HD 140283 (note that, 



-7- 



for comparison, the departure coefficients shown here is calculated with the the same stellar 
parameters as K94, while the final abundances of HD 140283 are determined with the stellar 
parameters given in Table [2]). It can be seen in Figured that our results confirm the conclu- 
sions of K94 that the departures from LTE for the resonance lines are due to a combination 
of overionization and optical pumping effects. However, the amplitude of the departures 
from LTE for our results (see Figured) are smaller than that of K94 (c.f. Figure 4 of K94). 
As described in § 14. 1[ the main differences in the boron atomic model between K94 and this 
work lie in: 

• excitation energies for some of the higher levels; 

• number of radiative b-b transitions considered; 

• photoionization cross- sect ions for the 4 P and 2 F° levels; 

• electron impact excitation cross-sections; 

• cross-sections of ionization induced by collisions with neutral hydrogen. 

The marginal differences in excitation energies and photoionization cross-sections of the 
higher levels cannot produce such large differences in the NLTE effects. The effects of the 
number of radiative b-b transitions considered was tested by excluding the transitions with 
/ < 0.001 as K94 did, but calculated results showed that there were no obvious differences. 
The differences in NLTE effects caused by the differences in electron impact excitation cross- 
sections was investigated in an extreme way, i.e., calculation was performed without the 
electron impact processes; in this case, the departures from LTE for our results increase a 
bit, but still are smaller than that of K94. Lastly, we checked the influences of the hydrogen 
collisional ionization on the NLTE effects of B I with an elaborate routine, which enabled 
us to treat the hydrogen collisional excitation and ionization separately. We first fixed 
the strength of hydrogen collisional excitation, and then calculated the NLTE effects with 
two different strengths of hydrogen collisional ionization (Sh = and 1, respectively). As a 
result, only marginal differences in the NLTE effects were found, which means that hydrogen 
collisional ionization process plays an insignificant role in the B I line formation. Similar 
conclusion was also derived by K94, who found the hydrogen collisional ionization rates to be 
several orders of magnitude smaller than the corresponding rates due to electron collisions. 
In a word, it's unlikely that the large differences in the NLTE effects between K94 and this 
work are due to the differences in atomic model. 

As pointed by K94, both the overionization and optical pumping effects are caused 
by the hot, non-local, ultraviolet radiation fields in the line formation regions, while the 
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background opacities in the NLTE calculation code used by K94 were just the continuous 
ones, i.e., the background line opacities were not included. To compensate this defect, 
K94 approximated the line blanketing effects on the photoionization processes by treating 
the transitions at fixed rates which were calculated from a specified radiation field. These 
radiation fields were produced directly by the adopted model atmosphere code, which had 
included the line blanketing via opacity sampling. However, such approximation was not 
applied to the radiative b-b transitions (except for some test runs), therefore the efficiency 
of the pumping processes would be exaggerated. In their subsequent work, KC96 used a 
revised version of NLTE calculation code which allows the inclusion of spectral lines in the 
background opacities. As a result, the departures from LTE decreased. For example, the 
NLTE abundance correction (A loge(B)@ = log e(B) NLTE -log e(B) LTE ) for the B I 2496.771 A 
resonance line decreased from 0.59 to 0.49 dex for HD 140283 after including background line 
opacities. However, the departures from LTE for the KC96's results are still larger than 
ours (as will be presented below, our NLTE B abundance correction for HD 140283 amounts 
only to 0.2 dex when neglecting collisions with neutral hydrogen). The NLTE calculation 
code used in this work has included m ost of the newe st continuous and line opacities, while 



some of them (the Fe I b-f opacities of iBautistal 119971 and the VALD atomic line data being 



the most prominent) were not available for KC96 more than one decade ago. It might be 
possible that KC96 used smaller background opacities than that adopted in this work and 
thus derived larger NLTE effects for B I. 

Another difference between the results of K94 and this work worthwhile to be noted 
is that K94 concluded that the departures from LTE would increase if including the b-b 
inelastic collisions with neutral hydrogen. This is completely contrary to our results that the 
line formation of B I will get closer to LTE if the strength of collisions with neutral hydrogen 
increases, as can be seen clearly in Figure [2j As the role of the hydrogen collisional ionization 
has been shown to be insignificant, the obvious changes in the NLTE effects should be mostly 
due to the hydrogen collisional excitation. 

Figure [3] shows the NLTE abundance corrections for the B I 2496.771 A resonance line as 
a function of stellar parameters. In general, the NLTE effects are large for hot, metal-poor, 
and low-gravity stars. This is in agreement with the results of KC96. 



5 loge(X) = log[iV(X)/7V(H)] + 12 
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Abundances and uncertainties 



5.1. Boron 

B abundances were derived by spectrum synthesis of the B I 2497 A resonance doublet 
using the IDL/FORTR AN SIU softwar e packa ge developed by Dr. J. Reetz. The input line 



list was extracted from iDuncan et al.l (ll998al ) with minor changes. We took the weaker 



2496.771 A component as the primary abundance indicator because the stronger 2497.725 A 
line is severely blended. Though the weaker line is also blended with a Co I line at 2496.708 A, 



(1997 


) and 


Primas et al. 


(1999) 



for the rest four stars, as reported by iGarcia Lopez et al.l (119981 ) . a combination of high 
effective temperature (T eff > 6000 K) and low metallicity ([Fe/H] < —2) prevented us from 
obtaining accurate B abundances with the B I 2496.771 A line, thus only upper limits were 
given. Figure H] shows the observed and synthetic spectra of the B I 2497 A resonance doublet 
for four program stars. 

The uncertainties associated with B ab undances were estimated following the same 
procedure described by IDuncan et al.l ( 119971 ) . Individual errors caused by uncertainties in 
effective temperature, surface gravity, metallicity, microturbulent velocity, photon statistics 
in observed spectrum, and continuum placement (a typical value of ±2% was adopted) were 
added in quadrature. The B abundances and final errors are given in Table [2j 

As mentioned abo ve, B abundances of our sample stars have be e n studied previousl y 
by other investigators ( IDuncan et al.l Il997t IGarcia Lopez et al.l Il998l ; iPrimas et al.lll999l ). 
Figure [5] shows the comparison of LTE B abundances between the literatures and this work for 
the 12 stars with determined B abundances. It can be seen that the agreement is reasonable 
for most of the stars. The only exceptio n is H P 184499, for which our B abundance is 
0.58 dex higher than that of IDuncan et al.l (Il997l). Such a l arge d ifference is mostly due to 
the different stellar parameters adopted by IDuncan et al. ( 1997 ) and by this work (their 
effective temperature and metallicity are 12 K and 0.4 dex lower than ours, respectively). 
If we adopted the same stellar parameters as IDuncan et al.l (119971 ) , then the difference in B 
abundance will decrease from 0.58 to 0.28 dex, which is comparable to the uncertainties. 



5.2. Oxygen 



As it is believed that boron production is related to oxygen directly, O abundance 
would be a preferred alternative to Fe in investigating the origin and evolution of B in the 
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Galaxy. There are several indicators for O abundance: the UV OH lines, the [O I] 6300 and 
6363 A forbidden lines, the O I 7774 A triplet, and the infrared (IR) vibration -rotation OH 



l ines, and different indicators may give different O abundances. As discussed in lNissen et al. 



(120021 ) and references therein, O abundances from the OH and O I lines are very sensitive 
to the adopted stellar parameters, such as the effective temperature and surface gravity; 
moreover, line formations in cool metal-poor stars are probably far from LTE for either the 
UV OH lines or the O I triplet. In contrast, [O I] forbidden lines are formed very close to 
LTE and O abundances from these lines are not sensitive to the adopted stellar parameters. 
Therefore, it is normally believed that [O I] forbidden lines are the most reliable indicators 
for O abundances, but the difficulty is that these lines are very weak in metal-poor dwarf 
and subgiant stars. 

For our sample stars, in the absence of reliable data for the [O I] forbidden lines, we used 
the O I triplet to determine their O abundances. Based on the equivalent widths of the O I 
7774 A triplet collected from the literatures, O a bundances were first calculated with log gf = 



0.369, 0.223, and 0.002 from Wiese et all f 



996) in LTE. Then NLTE corrections were applied 



according to the re sults of IFabbian et al.l ( 120091 ) which were calculat ed with new electron 
collisional data from lBarkleml ( 120071 ). We were aware that th e results of IFabbian et al.l (120091 ) 
were based on the MARCS models (IGustafsson et al.l 120081 ) . while our LTE O abundances 



were d etermined with the MAFAGS-OS9 models. Howe ver, we also noted that IFabbian et al. 



( 120091 ) found that, when using the ATLAS9 models of ICastelli fc Kuruczl ( 120031 ). resulting 
NLTE abundance corrections are very similar to that obtained using MARCS models for 
intermediately low metallicity ([Fe/H] > —2.5) stars. This is not unexpected because of 
the the similarity between ATLAS9 and MARCS models. As the MAFAGS-OS9 model 
used in this work is also quite sim ilar to the MARCS m odel atmosphere^], we suggest that 
the NLTE corrections de rived by IFabbian et al.l (120091 ) can be applied to our results too. 
The NLTE corrections of IFabbian et al.l (120091 ) were calculated with two different hydrogen 
collis ional strengths, nam e ly Sw = and 1. The latter was adopted in this work because 
both lAUende Prieto et al.l ( 120041 ) and iPereira et al.l ( 120091 ) found that the Su = 1 case could 
reproduce the solar center-to-limb variations better than S-r = 0. Moreover, as stated by 
Fabbian et al.l (120091 ). in the case of Sh = 1, derived [O/Fe] trend becomes almost flat below 
[Fe/H] ~ — 1, which is in better agreement with the results obtained form other O abundance 
indicators; and if Sh = is adopted, then [O/Fe] will decrease towards lower [Fe/H], which 
would open new questions. The final NLTE-corrected O abundances are given in Table [2] 



6 We have made comparisons between MAFAGS-OS9 and MARCS for a set of model atmospheres across 
the parameter space investigated here, and results show that the temperature, electron density, and gas 
pressure stratifications of the MAFAGS-OS9 models are very similar to that of the MARCS models. 
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(the reference solar O abundance is loge(O) = 8.77 from lTan et al.ll2009l ; we emphasize that 
the slopes of the Be(B)-0 abundance relationships, which are of our primary concern, are 
independent of the choice of the reference solar O abundance). 

The uncertainties of O abundances were estimated from the errors of the stellar pa- 
rameters and the equivalent widths. Typical errors of ±80 K in effective temperature and 
±0.15 dex in surface gravity both lead to an error of ±0.05 dex in O abundance. The uncer- 
tainties in iron abundance and microturbulent velocity nearly have no effect on O abundance. 
A typical error of ±3mA in equivalent width corresponds to an error of ±0.04 dex. In total, 
the typical error of O abundance is around ±0.08 dex. 



5.3. Beryllium 



Beryllium abundance can be a useful tool to understand the origin and evolution of 
boron. It is suggested that Be is only produced by spallation reactions, while B may be 
additionally yielded by the ^-process, therefore, the B/Be ratio can put some constraints on 
the production history of B. We computed the Be abundances from the equ ivalent widths 
of the B e II 3131.066 A resonance line (logqf = —0.468 from Kuruczl Il992l was adopted) 
given by iBoesgaard fc King (119931 ) and iBoesgaard et al.l ( 119991 ) for 13 of our sample stars. 
For the remaining three stars, Be abundances were adopted from the literatures but scaled 
to our stellar parameters. Beryllium abundances of all the stars were determined in LTE 
as NLTE effects for the B e II resonance lines are insignificant according to the results of 
Garcia Lopez et al.l ( 119951 ). 



The uncertainties of Be abundances mainly come from the errors of the stellar parame- 
ters and equivalent widths. The uncer t aintie s asso ciated with the e q uivale nt widths, which 
were not given by IBoesgaard fc King! ( 119931 ) and IBoesgaard et al.l (119991) . depend mostly 
on the placement of continuum. As discussed by iDuncan et al.l ( 119971 ) for boron, the con- 
tinuum is easier to define for more metal-poor stars, but the B lines are weaker; for stars 
with higher metallicity, the errors in continuum placement are larger but the B lines are 
stronger, so the combined effects are relatively consta nt over the metallic i ty ran ge investi- 
gated. This situation holds also for Be. We noted that IBoesgaard fc King! ( 119931 ) claimed a 
typical relative e r ror of about 1% (in the worst cases 2%) for their continuum placement, and 
Boesgaard et al.l ( 119991 ) estimated the typical Be abundance errors due to the placement of 
the continuum to be 0.02-0.06 dex. In this work, we adopted a conservative value of 0.1 dex 
as the typical uncertainty caused by the errors in continuum location. The Be abundances 
and final errors are listed in Table [2j 
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Discussion 



In this section we will investigate the implications for the origin and evolution of boron 
from the abundances of our sample stars. Before that it is necessary to check whether some 
of the stars are depleted in B because B can be destroyed in stars by fusion reactions at 
a relatively low temperature (about 5 x 10 6 K). This can be done by investigating their 
Li and Be abundances because the destruction temperature for Li, Be, and B increases in 
order. If Li or Be is not depleted, then B should not be depleted either. Among the 12 
stars with determined B abundances, five are ve ry lik ely depleted in Li (f ar below the Spite 
plateau) according to iGarcia Lopez et al.l ( 119981 ) and iPrimas et al.l ( 119991 ) . In the following 
some comments are given for these Li-depleted stars: 



HD 64 090 & HD 1844 99: though these two stars are obviously depleted in Li, their 
Be and B abundances match well with the Be-Fe(O) and B-Fe(O) trends (abundance 
trends in logarithmic plane, this meaning holds true for all the following denotations in 
the same style if not specified) determined by the Li-normal (thus Be- and B-normal) 
stars as can be seen in Figure [6j We conclude that their Be and B abundances are not 
depleted and will take them into account when investigating the origin and evolution 
of B. 

HD 106516 & HD 221377: both of them have only upper limits for their Li and Be 
abundances . Their B abundan ces are slightly lower than that of stars with similar 
metallicity. IPrimas et al.l (119991 ) also took these two stars as possible B-depleted stars. 
We will exclude them in the following discussion. 

BD +23° 3130: this is the most metal-poor star with determined B abundance in our 
sample. It is possibly a giant star with T c g = 5255 K and \ogg = 2.89. Its Be and B 
abundances are 0.37 and 0.25 dex lower, respectively, than that of HD 140283, which 
has very similar metallicity. It's hard to say whether this star is slightly depleted in Be 
or B, but for reliability, it will not be considered either in the following investigation. 



Another star worthy of being noted is HD 160617, which was concluded as a Li-normal 
but B-depleted s t ar (w hether it is Be-depleted depends on the adopted stellar parameters) 
by IPrimas et al.l (119981 ) . Both the Li an d Be abundance s of this star determined with our 



stellar parameters seem to be normal (seelTan et a 



work is 0.24 dex higher than that of IPrimas et al. 



2009). Its B abundance obtained by this 
jl998h which was determined with very 
similar parameters, but considering the uncertainties, they are still in reasonable agreement. 
However, we do note that Be and B abundances of HD 160617 are much smaller than that of 
HD 64090 which has very similar Fe abundance. This large differences lead to the relatively 
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large scatter at [Fe/H] ~ —1.8 on the Be-Fe and B-Fe trends as can be seen in Figure 
But on the Be-0 and B-0 diagram, the scatter is much smaller as shown in Figure [6b. This 
can be a direct evidence that Be and B production are directly related to O rather than 
Fe. Therefore, we take HD 160617 as a B-normal star and will include it in the following 
investigation on the B production and evolution. Finally, we have nine stars with normal Be 
and B abundances and the following discussion will be only based on these stars. 



6.1. Evolution of Be and B with metal abundances 



Figure E] shows the Be and B abundances against [Fe/H] and [O/H] for the nine Be- 
and B-normal stars. It is obvious that both loge(Be) and loge(B) increase approximately 
linearly with increasing [Fe/H] or [O/H]. Slopes and intercepts for the linear leat-squares fits 
to the abundance data are given in Table [3j for comparison, results from previous studies 
are also given. It can be seen that the slope of the Be-Fe trend obtai ned by this work agrees 
well with the result of iTan et al.l (120091 ). and is consistent with that of lSmilianic et al. (2009) 



within error bars. But our B- 



e tre nd is marginally steeper than that of IBoesgaard et al. 



( 119991 ) and lRich fe Boesgaardl (120091 ) . which is mainly due to the different metallicity ranges 



investigated. When the fitting is onl y restricted to the simila r met allicity range to thi s work , 
the slopes of the Be-Fe trends from IBoesgaard et al.l ( 119991 ) and iRich fe Boesgaardl (120091 ) 
are 1.17 and 1.07, respectively, which are consistent with our resul t. For the B-Fe trend, 



the slope obtained by this work is steeper than that of iDuncan et al.l ( 119971 ). This difference 
can be att ributed to the two lo west metallicity stars (BD+3°740 and BD — 13°3442) in the 
sample of IDuncan et al.l (J1997J), for which only upper limits of B abundances are derived 
by this work, and thus are not incl uded in our investiga tion of the B-Fe(O) relationships. 
It can be seen clearly in Figure 3 of IDuncan et al.l (119971 ) that, compared with the general 
B-Fe trend, these two stars show obvious enhancement of B abundances, and if they are 
excluded from the fitting, then the slope of the B-Fe trend will increase to 1.04, which is 
consistent with our result within uncertainties. For the Be-0 and B-0 trends, comparisons 
become more complicated because different work may use different O abundance indicators 
and/or different descriptions of the line formation (LTE versu s NLTE). The sl ope of the 
Be-0 trend obtained by this work agrees well with the result of ITan et al.l ( 120091 ). of which 
the O abunda nces were also d e termi ned using the O I triplet and had been corrected for 
NLTE effects. IBoesgaard et al.l ( 119991 ) adopted the averag e value s of O abundances based on 
the UV OH lines and the O I triplet; IRich fe Boesgaardl (120091 ) determined O abundances 
from the UV OH lines, therefore, a dire ct comparison of their results with that of this work 
is inappropriate. ISmiljanic et al.l ( 120091 ) used [a/H] as a proxy for [O/H] in the absence of 
O abundances for all their sample stars, but interestingly, the Be-a trend obtained by them 
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agrees well with the Be-0 trend of this work. Oxygen abundances adopted by lDuncan et al. 
( 119971 ) for most of their sample stars were based on the O I triplet and had been corrected 
for NLTE effects, but for the same reason as we have mentioned in the comparison of the 
B-Fe trend, their B-0 trend is flatter than that of this work. 

Both the slopes of our Be-Fe and Be-0 trends are closer to unity, which confirms the 
previous conclusions that Be should be mainly produced by primary process in the early 
Galaxy. As for B, the derived B-Fe and B-0 relationships depend somehow on the adopted 
strength of collisions with neutral hydrogen as we have discussed in § 14.21 As shown in 
Figure [3b, on one side, NLTE abundance corrections will increase if the strength of collisions 
with neutral hydrogen decreases; on the other side, NLTE abundance corrections increases 
with decreasing metallicity, therefore, the B L TE-Fe(0) and B NLTE SH=0 -Fe(O) trends shown 
in Figure [6] represent the two extreme cases, and the real situation should lie between them. 
In other words, the slopes of the B-Fe and B-0 relationships should be restricted to the 
intervals [1.07, 1.17] and [1.34, 1.46], respectively. In any case, the slopes of the B-Fe and B- 
O trends are closer to 1, which indicates that primary process is dominant in B production 
in the early Galaxy. In addition, there seems to be slight changes in slopes for both the 
Be-Fe(O) and B-Fe( O) trends around [F e /H] ~ —1 and [O/H] ~ —1, respectively, which 
was also reported by iGarcfa Lopez et all (119981 ). but there are too few stars in our sample 
to confirm it. 

It is believed that both Be and B can be produced by spallation reactions while B may 
be additionally produced by the //-process, therefore, B should increase more rapidly than 
Be (or at a similar speed as Be if th e cont ribution of the ^-process is neglectable). However, 
as first pointed by iBoesgaard et all (120041 ). after applying the NLTE abundance corrections 
of KC96, the B abundance in creases more slowly for halo st ars than Be does (see T ables 5 
and 6 of iDuncan et all Il997l and Figures 11, 13, and 14 of IBoesgaard et al] 120041 ). which 
is in contradiction with the theoretical models. For our results, as shown in Figure [6J this 
problem no longer exists no matter the strength of collisions with neutral hydrogen may be. 
Even in the extreme case (no collisions with neutral hydrogen), B increases with Fe or O 
at a very similar speed as Be does, and if the strength of collisions with neutral hydrogen 
increases, then B will increase with Fe or O a bit faster than Be. To further investigate 
whether our new result (B increase with Fe or O in a similar fashion to Be) are due to our 
smaller NLTE abundance corrections or to the LTE B abundances, we applied the NLTE 
abundance corrections of KC96 to our LTE B abundances. As a result, the slopes of the 
B-Fe and B-0 trends decrease to 0.97 and 1.22, respectively, which are lower than the slopes 
of the Be-Fe and Be-0 trends. In addition, a reverse test was also p e rform ed, i.e., we applied 
our NLTE corrections to the LTE B abundances of IDuncan et al.1 (ll997f ). and in this case, 
the slope of the B-Fe trend would be 0.86 (0.70 when the NLTE corrections of KC96 are 
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applied), which is still similar to the slope of the Be-Fe trend (0.85). Therefore, our result 
that B increases with Fe or O at a very similar speed as Be is mostly due to the smaller 
NLTE abundance corrections. 

Though good linear relationships between Be and Fe as well as between B and Fe 
are found for the metallicity range (—2.5 < [Fe/H] < —0.5 ) investigated here, it is not 
clear whether such trends will continue to lower metallicities. iPrimas et al.l ( 120001 ) reported 
the detection of Be in the very metal-poor star G 64-12 ([Fe/H] = —3.3) and found that 
Be abundance of this star is significantly higher than what expected from the general Be- 
Fe trend. This led them to suggest a possible flattening of Be abundances at very low 
metallicities. We also no te that, for BD+3°740 an d BD — 13°3442, only up per limits of B 



abundances a re found by iGarcfa Lopez et al.l (119981 ) and by this work, while iDuncan et al. 



( 119971 . Il998bf ) claimed the detection of B in these two stars. And inter estingly, B abundances 
of BD +3°740 and BD -13°3442 given bv buncan et al.l (1l997i . ll998bl ) also seem to be higher 
than what expected from the gener al B-Fe trend (both of these two stars are at the metallicity 
end of their sample). Moreover, lAsplund et al.l (120061 ) reported the detection of 6 Li (at 
the significance level > 2a) in nine metal-poor stars. Their most noteworthy result is the 
detection of an abnormally high 6 Li abundance in the very metal-poor star LP 815-43, which 
cannot be explained either by the standard BBN or by the Galactic cosmic-ray spallation 
and a-fusion reactions. One common solution to all the three problems described above can 
be non-standard BBN, which could produce much more primordial light elements than the 
standard BBN. However, if ignoring the problem of high 6 Li abundance at very low metallicity 
(whether 6 Li is really d etected in some metal-poor stars is still disputable, for example, 
Garcia Perez et al.ll2009l aimed to determine 6 Li/ 7 Li for a sample of metal-poor stars, and 
only found that the 6 Li/ 7 Li ratios are comparable to or even lower than the associated 
uncertainties), there may be no need to invoke the non-standard models o f BBN to explain 
the possible Be and B "plateau" as pointed by IVangioni-Flam et al.l (119981 ) . They proposed 
that in the very early Galaxy, if supernovae of all masses paly roles in the production of Be(B), 
then Be(B) should continue to decrease with decreasing [Fe/H]; otherwise, if only the most 
massive stars (initial mass larger than 60 M ) contribute to the Be(B) yields, then Be(B) 
should be roughly constant. They suggested that the predicted di fference in the behavior o f 
Be(B) could be tested by observations at [Fe/H] < —3. Recently, iRich &: Boesgaardl (120091 ) 
investigated the Be abundances of 24 very metal-poor stars with [Fe/H] from —2.3 to —3.5 
and found changes in slopes for the Be-Fe trend at [Fe/H] ~ —2.2 as well as for the Be-0 
trend at [O/H] ~ —1.6. In general, their results show that Be increases with Fe(O) much 
more slowly below [Fe/H] ~ —2.2 ([O/H] ~ —1.6) than above that. Though there is no 
indication of a Be plateau, the four most metal-deficient stars ([Fe/H] < —3) do show a Be 
enhancement with respect to Fe at the la level, which is qualitatively consistent with the 
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predictions of the light-element product ion model involved with only the most massive stars 
proposed by lVangioni-Flam et al.l ( 119981 ). To clarify this problem, more accurate data on Be 
and B abundances at very low metallicities are needed. 



6.2. The B/Be ratio 



Figure [7] shows the B/Be ratio as a function of metallicity for the nine Be and B- normal 
stars. Linear least-squares fits to the data give the following relationships: 

log(B LTE /Be) = (0.10 ± 0.14) [Fe/H] + (1.25 ±0.21) 
log(B NLTEiSH=0 /Be) = (0.01 ±0.12) [Fe/H] ± (1.22 ±0.19) 

It is clear that, in LTE case, B/Be increases slightly with increasing metallicity, while in the 
case of NLTE but without neutral hydrogen collisions, B/Be is almost constant. Again, the 
real situation for the B/Be-Fe trend should lie between the two extreme cases, but considering 
the uncertainties, the B/Be ratio is roughly constant over the metallicity range investigated 
here. The average B/Be ratio for the nine stars should fall in the interval [13 ± 4, 17 ± 4], 
of which the lower and upper limits correspond to the two extreme cases re spectively This 
is in reasonable agreemen t with the results 15 ± 3 of iDuncan et al.l ( 119971 ) and 20 ± 10 of 
Garcia Lopez et al.l (119981 ) . Our average B/Be ratio is consistent with the predicted value of 
spallation process, which depends on the compositions and energy spectra of c osmic rays (for 



exam ple, a typical spallation-produced B/Be ratio of about 14 is derived by iRamaty et al. 
2000l ). However, if assuming that B is purely produced by spallation reactions and the 



11 B/ 10 B ratios in metal-poor stars are the same as the mete oric value, then a minimum 
B/Be ratio of about 20 is required (jVangioni-Flam et al.lll996l ). This value is slightly larger 



than o ur result, which might indicate that B is also produced by the //-process. IRamaty et al. 



( ]2000l ) showed that even though the production histories of the spallation-produced B and 
Be and the neutrino-produced B are different, B/Be can still remain essentially constant as 
a function of [Fe/H]. However, the 11 B/ 10 B ratios in metal-poor stars need not to be the 
same as the meteoric value, especially considering the fact that our average B/Be ratio in 
metal-poor stars is lower than the solar meteoric value of about 23. Anyway, an accurate 
measurement of the 11 B/ 10 B ratios in metal-poor stars is strongly desirable. 



7. Summary 



In this work, we have investigated the NLTE line formation of neutral boron in the 
atmospheres of cool stars. Our results confirm the conclusions of K94 and KC96 that NLTE 



-17- 



effects for the B I resonance lines, which are due to a combination of overionization and 
optical pumping effects, are significant for hot, metal-poor, and low-gravity stars. However, 
the amplitude of the departures from LTE found by the present study are smaller than that 
of KC96, which might be because that KC96 used smaller background opacities than that 
adopted in this work. In addition, our calculation shows that departures from LTE for the 
B I resonance lines decrease with increasing strength of collisions with neutral hydrogen, 
which is completely contrary to the result of K94. 

The NLTE line formation results for B I have been applied to the determination of B 
abundances for a sample of 16 metal-poor stars, for which the Be and O abundances are 
also determined. No matter the strength of collisions with neutral hydrogen may be, both 
the slopes of the Be-0 and B-0 trends determined by the nine Be- and B-normal stars are 
closer to unity, which confirms the previous conclusions that primary process is dominant in 
the production of Be and B in the early Galaxy. Most importantly, our results show that, B 
increases with Fe or O at a very similar speed as, or a bit faster than Be does, which is in 
accord with the theoretical models. The B/Be abundance ratios are roughly constant over 
the metallicity range —2.5 < [Fe/H] < —0.5. The average B/Be ratio falls in the interval 
[13 ± 4, 17 ± 4], which is consistent with the predictions of spallation process. However, if 
assuming the isotopic 11 B/ 10 B ratios in metal-poor stars to be the same as the meteoric 
value, then a minimum B/Be ratio of about 20 is required by the pure spallation process, 
therefore the //-process may be needed to reconcile this problem. But whether the 11 B/ 10 B 
ratios in metal-poor stars are the same as the meteoric value is still unclear, thus an accurate 
measurement is strongly desirable, though rather challenging. 
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Fig. 1. — Grotrian diagram for B I (the ground state of B II is also shown). The continuous 
lines denote the radiative bound-bound transitions considered in this work. 
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Fig. 2. — Departure coefficients for the populations of first eight energy levels of B I and the 
ground state of B II for the representative metal-poor star HD 140283 calculated with the 
stellar parameters from K94. 



-24 - 




Fig. 3. — Dependence of NLTE abundance corrections for the B I 2496.771 A resonance 
line on (a) effective temperature, (b) metallicity, and (c) surface gravity. Squares, triangles, 
and diamonds represent results calculated with different strength of collisions with neutral 
hydrogen, i.e., Sh = 0, 0.1, and 1 respectively. 
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Fig. 4. — Spectrum synthesis of the B I 2497 A resonance doublet for four program stars. The 
observed spectra are represented by photon statistics error bars. The solid line is the best- 
fitting synthesis, and the dotted lines are synthetic spectra with B abundances of ±0.2 dex 
relative to the best synthesis. 
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Fig. 5. — Comparison of B abundances (LTE) between the literatures and t his work 
for the 12 stars with dete rmined B abundances: d iamon ds, iDuncan et al.l (119971 ): square, 
Garcia Lopez et al.l ( 119981 ); triangles, iPrimas et al.l ( 119991 ). 
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Fig. 6. — Be and B abundances against (a) [Fe/H] and (b) [0/H] for the nine Be- and 
B-normal stars. Filled triangles are Be abundances, filled circles are LTE B abundances, 
and filled squares are NLTE B abundances calculated without neutral hydrogen collisions 
(Sh = 0). Solid, dotted, and dashed lines are the best linear fits to the Be-Fe(O), Blte-Fc(O), 
and B NLTEi s H= o-Fe(0) relationships respectively. 
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Fig. 7. — B-to-Be abundance ratio as a function of [Fe/H] for the nine Be- and B- normal 
stars (filled circles: B abundances calculated in LTE; filled squares: B abundances calculated 
in NLTE without neutral hydrogen collisions, i.e., Sn = 0). Dotted and solid lines represent 
the best linear fits to the BLTE/Be-Fe and BNLTE,s H =o/Be-Fe relationships, respectively. 
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Table 1. Stellar parameters 



Stf fir 


J- eff 


iu & y 


\Fp /Hi 
[re/nj 


P 
s 


Ref 




(K) 






( \r m c ~ \ 
I 1V111 J 




HD 19445 


6009 


4.40 


-2.02 


1.6 


1, 2 a 


HD 64090 


5525 


4.57 


-1.80 


2.0 


5 


HD 76932 


5890 


4.12 


-0.89 


1.2 


4 


HD 84937 


6356 


4.00 


-2.16 


1.7 


1, 3 a 


HD 94028 


5925 


4.19 


-1.51 


1.5 


2 


HD 106516 


6115 


4.28 


-0.80 


1.7 


5 


HD 140283 


5725 


3.68 


-2.41 


1.5 


4 


HD 160617 


5940 


3.80 


-1.78 


1.5 


4 


HD 184499 


5820 


4.07 


-0.59 


1.2 


5 


HD 194598 


5980 


4.27 


-1.16 


1.6 


2 


HD 201891 


5900 


4.22 


-1.07 


1.2 


2 


HD 221377 


6330 


3.81 


-0.78 


2.0 


5 


BD +3°740 


6290 


4.02 


-2.55 


1.5 


2, 3 a 


BD +23°3130 


5255 


2.89 


-2.42 


1.7 


6 


BD +26°3578 


6280 


3.93 


-2.25 


2.0 


2 


BD -13°3442 


6390 


3.88 


-2.66 


1.4 


3 



References. 



11 lKornetal.1 fl2003h: (2) 



Gehren et all (bood) : (3) iMashonkina et al.l (12008^ : (4) 



Tan et all (120091 ); (5) determined with ELODIE spectra 



by this work; (6) determined with UVES spectra by this 
work. 



a The average value adopted. 
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Table 2. Abundance results 



Star [O/H] Ref. loge(Be) Ref. <r(Be) loge(B) loge(B) NLTE a(B) 

NLTE LTE LTE S H = 5 H = 0.1 S H = 1 



HD 19445 


-1 


40 


8 


-0 


50 


7 





13 







56 







74 







67 







59 





19 


HD 64090 


-1 


08 


2 


-0 


02 


7 





15 




1 


24 




1 


33 




1 


27 




1 


25 





13 


HD 76932 


-0 


43 


11 





88 


7 





14 




2 


03 




2 


08 




2 


07 




2 


05 





14 


HD 84937 


-1 


69 


2 


-0 


85 


7 





12 


< 





43 


< 





73 


< 





72 


< 





55 






HD 94028 


-0 


95 


4 





33 


7 





13 




1 


24 




1 


35 




1 


31 




1 


26 





16 


HD 106516 


-0 


30 


1 


< -0 


58 


5 








1 


38 




1 


45 




1 


43 




1 


41 





15 


HD 140283 


-1 


77 


8 


-0 


96 


7 





14 







07 







27 







18 







11 





22 


HD 160617 


-1 


38 


8 


-0 


41 


12 





12 







51 







66 







62 







55 





20 


HD 184499 


-0 


27 


3 


1 


02 


7 





16 




2 


20 




2 


23 




2 


23 




2 


22 





17 


HD 194598 


-0 


77 


11 





18 


7 





14 




1 


38 




1 


46 




1 


44 




1 


40 





14 


HD 201891 


-0 


64 


9 





51 


7 





14 




1 


72 




1 


79 




1 


76 




1 


73 





15 


HD 221377 


-0 


63 


1 


< -0 


92 


3 








1 


41 




1 


50 




1 


50 




1 


48 





19 


BD+3°740 


-2 


16 


2 


-1 


14 


7 





13 


< 





44 


< 





77 


< 





74 


< 





54 






BD+23°3130 


-1 


66 


4 


-1 


33 


6 





14 







18 







02 







12 







16 





25 


BD +26°3578 


-1 


73 


8 


-0 


77 


7 





13 


< 





24 


< 





52 


< 





52 


< 





36 






BD-13°3442 


-2 


22 


10 


-1 


29 


7 





12 


< 





23 


< 





55 


< 





52 


< 





33 
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Table 3. Results of linear least-squares fits to the Be and B abundance data 





X 


y 


Slope 


Intercept 






This work 




log 


;e(Be) 


[Fe/H] 


1.11 ±0.10 


1.74 ± 0.16 


log 


; e(B) LTE 


[Fe/H] 


1.17 ± 0.13 


2.96 ± 0.19 


log 


,'e(B)NLTE,£ 


iH =o [Fe/H] 


1.07 ±0.12 


2.92 ± 0.19 


log 


;e(Be) 


[0/H] 


1.33 ±0.13 


1.40 ± 0.14 


log 


;e(B) LTE 


[0/H] 


1.46 ±0.16 


2.63 ± 0.17 


loge(B) NLT E,s H =o [0/H] 


1.34 ±0.15 


2.62 ± 0.16 






Rich & Boessraard (2009) 




log 


;e(Be) 


[Fe/H] 


0.92 ±0.04 


1.41 ± 0.09 


log 


;e(Be) 


[O/H] 


1.21 ±0.08 


1.32 ±0.14 






Tan et al. 


(20091 




log 


;e(Be) 


[Fe/H] 


1.10 ±0.07 


1.63 ±0.10 


log 


;e(Be) 


[0/H] 


1.30 ±0.08 


1.19 ±0.08 






Smilianic et al. ( 20091 




log 


;e(Be) 


[Fe/H] 


1.24 ±0.07 


1.62 ±0.08 


log 


;e(Be) 


[«/H] 


1.36 ±0.08 


1.38 ±0.07 






Boeseaard et al. (1999) 




log 


;e(Be) 


[Fe/H] 


0.96 ±0.04 


1.41 ±0.03 


loge(Be) 


[0/H] 


1.45 ±0.04 


1.31 ±0.04 






Duncan et al. (19971 




log 


; e(B) LTE 


[Fe/H] 


0.96 ±0.07 


2.59 ±0.13 


log 


;e(B) LTE 


[0/H] 


1.21 ±0.13 


2.25 ±0.16 



